Observation of nonlinear dispersion relation and spatial 
statistics of wave turbulence on the surface of a fluid 
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We report experiments on gravity- capillary wave turbulence on the surface of a fluid. The wave 
amplitudes are measured simultaneously in time and space using an optical method. The full 
space-time power spectrum shows that the wave energy is localized on several branches in the wave- 
vector- frequency space. The number of branches depend on the power injected within the waves. The 
measurement of the nonlinear dispersion relation is found to be well described by a law suggesting 
that the energy transfer mechanisms involved in wave turbulence are not only restricted to purely 
resonant interaction between nonlinear waves. The power-law scaling of the spatial spectrum and 
the probability distribution of the wave amplitudes at a given wave number are also measured and 
compared to the theoretical predictions. 

PACS numbers: 47.35.-i,05.45.-a,47.52.+j 
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Wave turbulence concerns the study of the statistical 
and dynamical properties of a set of numerous nonlin- 
ear interacting waves. It is an ubiquitous phenomenon 
observed in various situations from spin waves in solids, 
internal or surface waves in oceanography up to plasma 
waves in astrophysics (for recent reviews see ll|,l2|). Wave 
turbulence theory, also called weak turbulence, predicts 
a wave energy cascade through the scales that can be 
derived analytically in nearly all fields of physics involv- 
ing weakly nonlinear interacting waves in infinite systems 
[3[. However, few well-controled laboratory experiments 
have been performed so far, and show partial agreement 
with the theory jl|, |2| . While most in situ or laboratory 
measurements involve time signals at a fixed location, 
theoretical predictions often concern the Fourier space. 
An important challenge is thus to get a space-time mea- 
surement of the turbulent wave amplitudes (as recently 
achieved for elastic wave turbulence [j] ) , and thus to have 
a better understanding of the elementary dynamical pro- 
cesses involved in the energy cascade. Concerning wave 
turbulence on a fluid, previous results involve either 2D 
spatial measurements but not resolved in time (oceanog- 
raphy [5| and laboratory experiments [8|) or resolved 
in time but restricted to ID space [6|, |7|. Here, we inves- 
tigate 2D spatial and temporal statistics of wave turbu- 
lence on the surface of a fluid by using an optical pro- 
filometry technique. We perform a Fourier analysis of 
movies of the free-surface deformation and focus notably 
on the nonlinear dispersion relation. 

The experimental setup consists of a tank (46 cm x 
36cm) filled with water (7cm deep). Surface waves are 
generated by the horizontal motion of two plunging rect- 
angular wave makers (19 cm in width and 2 cm in depth). 
They are located at two corners of the same longest side 
of the tank, the vibration directions being perpendicular 
to each other [ll|. The wave makers are driven by two 
electromagnetic shakers submitted to a random forcing 



within a narrow low-frequency band (typically from 1 to 
4 Hz). Typical maximal crest-to-trough wave amplitude 
ranges from 1 mm to 1.5 cm, and the wave mean steep- 
ness (ratio of crest-to-trough amplitude to its duration) 
ranges from 0.2 up to 3.3cm/s. This latter value corre- 
sponds to an injected power, P, 600 times greater than 
its value at the minimum forcing amplitude. This en- 
ables to access to linear, weakly and strongly nonlinear 
wave regimes. A Fourier transform profilometry method 
[Sl, llO| provides the temporal evolution of the vertical 
deformation of the free-surface of the fluid over a signif- 
icant spatial zone of the tank. Namely, a fringe pattern 
(wavelength A/ = 2.6 or 5.2 mm) is projected on the fluid 
surface by a video projector. When waves are generated, 
the vertical displacement of the free-surface leads to a 
phase shift of the pattern that is recorded by a camera. 
The deformation of the fluid surface r]{x, y^t) is then re- 
covered by a 2D phase demodulation of each image of 
the recorded movie [9|, llO|. Movies are recorded with 
1600 by 1200 pixels at facq = 50 or 60 Hz during roughly 
1 minute. The size of the recorded image is 25 x 19 cm^. 
To improve the contrast of the projected fringes on the 
fluid surface, a high concentrated white dye is added to 
the water bulk at an optimum concentration of 0.5% v/v 
[lO[ . The surface tension of this dyed water is measured 
to be 7 = 32 ±1 mN/m. Spatial and temporal resolutions 
of the measurement are 3 A/ and 2/ facq-, its linearity being 
insured for waves with sharp slopes up to 10 [10|. Pos- 
sible underneath hydrodynamic turbulence generated in 
the bulk by the wave makers does not play a significant 
role on wave turbulence (similar results are found when 
the immersed length of the wave maker is changed) . 

The vertical velocity of the fluid surface vix^y^t) is 
obtained by differentiating the wave height movie in 
time. The full space-time power spectrum of the veloc- 
ity £^(k, /) (a function of both the wave vector k and 
the frequency /) is then computed from multidimen- 
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FIG. 1: (Color online) Space-time spectrum E{k, f) of the 
vertical velocity of surface waves: moderate (a) and strong 
(b) injected powers [P^^^ = 5.2 (a) and 24.7 (b) in arb. units]. 
Forcing: 1-4 Hz. Colors are log scaled. Solid white lines are 
^N{k) with iV = 1, 2, and 3 (see text). Slope of the dotted 
line corresponds to a constant phase velocity uj{k)/k. 



sional Fourier transform. By integrating £^(k, /) over 
all directions of k, one obtains the velocity spectrum 
E{k = ||k||,/) displayed in Fig. [1] for moderate and 
strong forcings. We observe that the energy injected at 
low frequencies cascades through the scales and is mainly 
localized on several branches in the {uo = 27r/, k = 27r/A) 
space. At low forcing amplitude (not shown here), only 
one branch occurs that corresponds to the linear gravity- 
capillary relation dispersion a;(/c) = y^^+(7/p)^7 with 
g = 9.81 m/s^ the acceleration of gravity, p = 1000 
kg/m^ the fluid density. When the forcing is increased 
(see Fig. [TJi), a secondary branch appears below the lin- 
ear dispersion relation (LDR). This branch is found to 
be well described by flN{k) = ^/gNk + {'y/p)k^/N with 
N = 2 with no adjustable parameter (see solid line). At 
higher forcing (see Fig. [TJd), a third branch appears fol- 
lowing r^3(/c). Thus, as the power injected in the wave 
system increases, the nonlinear wave interactions redis- 
tribute the wave energy on N branches govern by QN{k)^ 
the nonlinear dispersion relation (NLDR). These sec- 
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FIG. 2: (Color online) Nonlinear dispersion relation k{f) com- 
puted from the lines of maximum energy oi E{k, f) for differ- 
ent forcings [P^/^ = 1 (+), 5.2, (o), 10.5 (x), 24.7 (*)]. Solid 
thick lines around each branch correspond to the branch width 
averaged for all forcings. Solid lines are Hat (A:) with A^ = 1, 
2, and 3 (see text - same as in Fig. [1]). Inset: snapshot of the 
wave amplitudes at strong forcing (P^^^ = 24.7). 



ondary branches arise from the propagation of harmonics 
(Nuj^Nk) superimposed on a carrier long wave (uj^k) and 
propagating with the phase velocity of the carrier (see be- 
low and Fig.fljD). Note also that ^N{k) = Nujj.^]^. Thus, 
at a fixed k^ corresponds N peaks (cj/c*, 2cj/e*/2 7 3cj/e^/3, 
• • • ) in a frequency Fourier spectrum, i.e. a horizontal 
slice of Fig. [H This is consistent with a two-peak fre- 
quency spectrum reported in a numerical simulation jl2| . 

At weak forcing, one observes linear gravity-capillary 
waves of gentle amplitudes that mix together. At strong 
forcing, steep long waves occur with sharp crest-ridges 
(see inset of Fig. [2j). Near the crests of these waves, high 
order harmonics are generated: small gravity-capillary 
waves superimposed on the long wave are observed (see 
also [l^). These harmonics are called bound waves since 
they do not propagate with their own phase velocity but 
with the one of the carrier long wave [13| , and thus leads 
to harmonics flN{K) of velocity Qn{K)/K = uj{k)/k 
where K = Nk. They thus do not obey the linear dis- 
persion relation which is consistent with the observation 
of secondary branches of the NLDR. This shows that 
other mechanisms than purely resonant wave interaction 
should be taken into account to describe the energy trans- 
fer across scales in wave turbulence. 

Let us now focus on the effect of the injected power P 
on the location and the width of branches QN{k). The 
maximum amplitude of each branch of the velocity spec- 
trum E{k^ /) is extracted using the maximum of a Gaus- 
sian fit with respect to /c at a fixed /. For different P, 
the lines of maximum energy of each branch are shown 
in Fig. [2] (symbols). Whatever the branch, the localized 
energy line is found to be independent of P: no measur- 
able shift of these lines occurs in the {k^u) space. This 
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FIG. 3: (Color online) Space-time spectrum E{kx, f) of veloc- 
ity at P^/^ = 10.5 located at ky = 0. Same forcing bandwidth 
as in Fig. [ll Inset: space spectrum E{kx,ky) located at / =5 
Hz (a) and 9 Hz (b). Dashed lines: forcing directions. Log 
scaled colors are different for each plot. 



differs from recent observation reported in elastic wave 
turbulence [J] or in simulation 12| . The widths of these 
branches are also plotted in Fig. [2l The width is de- 
fined by the rms value of the Gaussian fit. Whatever 
the branch, no significant evolution of the width is found 
when P is increased. The width is also independent of the 
branch number within our experimental accuracy. The 
typical width [A(A~^),A/] centered on a point (A~^,/) 
of the NLDR is roughly (8 m~^, 1.5 Hz), and could be 
ascribed to the typical nonlinear scale of wave mixing. 
However, one should be careful since A(A~^) is close to 
the resolution of the discrete Fourier Transform (~ in- 
verse of the image size ^ 5 m~^). To sum up, the wave 
energy is redistributed on different branches of the NLDR 
of width that is independent of P and of N. 

Figure [3] shows different views of the full space-time 
Fourier spectrum of the velocity £^(k, /). Main figure is 
E{kx^f)^ a slice at /c^y = of the spectrum along the 
X-axis. Figures [3^ and ^jp show £^(k) at two fixed fre- 
quencies (5 and 9 Hz). In Fig. [3^, the forcing reminisence 
appears as strong peaks in the direction of the wave mak- 
ers (see dashed lines). Due to wave turbulence, energy 
cascades across scales as can be seen in the two con- 
tinuous branches in the main Fig. [3l Simultaneously, a 
spread of energy across angles is responsible for the con- 
tinuous circles observed in the insets: the two concentric 
circles in Fig. [Sb correspond to two wave numbers given 
by both branches of the NLDR, whereas in Fig. [3^ one 
rather observes a single disc due to the overlapping of 
both branches of non-zero width (see above). Although 
not fully isotropic, this angular redistribution of energy is 
due to wave turbulence and is also visible as the symme- 
try between positive and negative kx in the main Fig. [3l 
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FIG. 4: (Color online) Spatial power spectrum E{k) of the 
velocity for P^/^ =1, 5.2, 10.5, 24.7 (from bottom to top). 
Dashed lines have slopes —4.6, —4.5 and —4.2 (from bottom 
to top). Dot-dashed line: A^^ = 25.2 m~^ (see Fig. 5). 



The space spectrum E{k = ||k||) of the velocity is then 
computed by summing the 3D space-time spectrum of 
£^(k, /) over all the directions of k and over /. Figured 
shows E{k) when the forcing is increased. At high enough 
forcing, E{k) is found to be scale invariant as expected 
for wave turbulence. The inert ial range increases with 
the forcing, and E{k) ~ k~^ with n c^ 4.2 over almost 
one decade in k corresponding to A ~ few cm. Note that 
n does not depend strongly on the forcing. Since one 
cannot compute E{f) from £^(k, /) in a wide range of / 
due to the strong steepness of the spectrum, one per- 
forms single point temporal measurements that shows 
a strong dependence of the power-law frequency expo- 
nent of the velocity spectrum on the forcing as already 
reported (typically from f~^ to /~^) [ll|, ll5|- These 
scalings suggest that the change of variable k ^ f us- 
ing the LDR to estimate E{k) from E{f) is not valid in 
temporal measurements in hydrodynamics wave turbu- 
lence when strong nonlinear waves are involved. Indeed, 
this would lead to an estimated velocity spectrum from 
k~^ to /c~^/^. Moreover, our velocity spectrum scaling, 
E{k) ~ /c~^"^, cannot be described by any of existing 
theories of wave turbulence taking into account either 
the presence of random-phased weakly nonlinear waves 
[Etheo{k) ^ P^/^/c~^/^] [16] or the dominance of coher- 
ent sharp wave crests [Etheo{k) ~ k~^ to k~^] [7[. It is 
known numerically that the spatial spectrum exponent 
can change in case of anisotropy [l7|. To confirm that 
our results do not depend on anisotropy, one computes 
E{k) by summing £^(k, /) over different {kx , ky) space 
regions: (i) over {kx ^ ky > 0) where isotropy is observed, 
(ii) over {kx , /c^ < 0) where anisotropy occurs due to the 



forcing (see insets of Fig. [3]), and (iii) all over {kx , ky) 
space. The spatial spectrum E{k) ~ k~^ computed over 
these different regions leads to z = 4.3 (i), 3.8 (ii) and 
4.2 (iii). Thus, the anisotropy does not play a significant 




of high Fourier mode amphtude is consistent with ID spa- 
tial measurements [7], simulations 12| and theory [18|. 
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FIG. 5: (Color online) PDF of the Fourier wave amplitude, 
|7}fcjVc^|r7fcj2, at k^ = 27t/X^ with A^^ = 25.2 ± 0.6 m"^ 
(see Fig. \^ for two forcings P^/^ = 5.2 (o) and 10.5 (*). 
(I^fc^l^) = 1-1 (o) and 1.3 (*). Solid lines have slopes —0.46 
and —0.58. Curves have been shifted vertically for clarity. 
Inset: \fik,\^ vs. time. A"^ = 25.2 ± 0.6 m-\ P^/^ ^ 24.7. 



In conclusion, we have reported 2D spatial statis- 
tics of wave turbulence on the surface of a fluid. The 
power spectrum, the nonlinear dispersion relation and 
the PDF of the Fourier modes show strong effects of non- 
hnear waves involved in wave turbulence. This suggests 
that energy transfer mechanisms are not only restricted 
to resonant interactions between nonlinear waves, but 
also involve the formation of localized nonlinear struc- 
tures (sharp-crested gravity waves) and of bound gravity- 
capillary waves. The wave spectrum scalings emphasize 
that the transition from /c-space to cj-space cannot be 
done according to the linear dispersion relation as usu- 
ally performed in wave turbulence experiments. 
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role on the estimation of the spatial spectrum exponent. 
Finally, the probability density function (PDF) of the 
wave amplitude rjix^y) is found to be roughly Gaussian 
whatever the forcing. One also computes the PDF of 
the Fourier amplitude \f}k^\^ of a wave component at a 
given wave number k-^. As shown in Fig. IH we choose 
k^ = 27r/A^ in the gravity regime with A~^ = 25.2 ± 0.6 
m~^ corresponding (using the LDR) to a frequency of 6.5 
Hz above the forcing ones. For each image, the value of 
\fjk^\'^ is extracted by averaging on 42 amplitudes found 
on a k-space ring of radius k^. Iterating for all images 
leads to the temporal evolution of the Fourier ampli- 
tude of the mode k^ as shown in Fig. [5l This signal 
is strongly erratic and bursts of random large-amplitude 
occurs. Similar random bursts of Fourier amplitude has 
been reported in simulations [12], these bursts being cor- 
related with phase jumps underlying strong nonlinear ef- 
fects [l2[. Although we are not able to measure the phase, 
this similarity is consistent with our above results under- 
lying strong nonlinear effect. The PDF of the Fourier 
amplitude \fjk^ P, rescaled to its rms value cr|^^ p, is then 
plotted in Fig. [5] for two forcings. At low forcing, the 
PDF is roughly exponential as expected for random and 
uncorrelated waves. At higher forcing, the PDF remains 
Gaussian up to three standard deviations, whereas its 
tail shows a slight departure from this Gaussian. Al- 
though more statistics are needed to characterize more 
deeply the PDF tail, this anomalously large probability 
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